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ABSTRACT 

The late stages of evolution of the primordial circumstellar disks surrounding young stars are poorly 
understood, yet vital to constrain theories of planet formation. We consider basic structural models 
for the disks around two ^lOMyr-old members of the nearby RCrA association, RX J1842. 9-3532 and 
RX J1852. 3-3700. We present new arcsecond-resolution maps of their 230 GHz continuum emission 
from the Submillimeter Array and unresolved CO(3-2) spectra from the Atacama Submillimeter Tele- 
scope Experiment. By combining these data with broadband fluxes from the literature and infrared 
fluxes and spectra from the catalog of the Formation and Evolution of Planetary Systems (FEPS) 
Legacy program on the Spitzer Space Telescope, we assemble a multiwavelength data set probing the 
gas and dust disks. Using the Monte Carlo radiative transfer code RADMC to model simultaneously the 
SED and millimeter continuum visibilities, we derive basic dust disk properties and identify an inner 
cavity of radius 16 AU in the disk around RX J1852. 3-3700. We also identify an optically thin 5 AU 
cavity in the disk around RX J1842. 9-3532, with a small amount of optically thick material close to 
the star. The molecular line observations suggest an intermediate disk inclination in RX J1842. 9-3532, 
consistent with the continuum emission. In combination with the dust models, the molecular data 
allow us to derive a lower CO content than expected, suggesting that the process of gas clearing is 
likely underway in both systems, perhaps simultaneously with planet formation. 
Subject headings: circumstellar matter — planetary systems; protoplanetary disks — stars: pre-main 
sequence — stars: individual (RX J1842.9-3532, RX J1852.3-3700) 



1. INTRODUCTION 

One of the key problems in planet formation is un- 
derstanding how the reservoir of planet-forming mate- 
rial - the disk of gas and dust around a young star - 
evolves with time. Perhaps the most compelling stage 
in the evolution of protoplanetary disks is the "tran- 
sitional" stage during which gas and dust are cleared 
from the system (e.g. Strom et al. 1989; Skrutskie et al. 
1990). This stage appears to be either rapid or rare, 
since fewer than 10% of low- to intermediate-mass stars 
in young star-forming regions are typically observed to 
be transitional systems (e.g. Silverstone et al. 2006; Cieza 
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et al. 2007; Uzpen et al. 2008), perhaps with some depen- 
dence on stellar mass (e.g. Currie & Sicilia-Aguilar 2010). 
Transition disks are traditionally identified empirically 
by a deficit of mid-infrared dust excess over the stel- 
lar photosphere relative to comparable young systems. 
This deficit is associated with a lack of hot dust close 
to the star. The advent of the Spitzer Space Telescope 
has revolutionized the quality and quantity of available 
data on the dust content of young stellar systems, partic- 
ularly transitional objects (see, e.g., Najita et al. 2007; 
Currie et al. 2009; MuzeroUe et al. 2010). One of its 
many contributions has been to reveal a new class of 
gapped or "pre-transitional" disks, in which an anoma- 
lously faint mid-infrared spectrum is bracketed by sub- 
stantial near- and far-infrared excesses (Espaillat et al. 
2007, 2008, 2010). Follow-up of S'piteer-identified transi- 
tional systems with high spatial resolution observations 
of continuum emission at millimeter wavelengths has led 
to the confirmation that mid-IR spectral deficits are as- 
sociated with a lack of long- wavelength emission from the 
disk center (Calvet et al. 2002, 2005; Hughes et al. 2007, 
2009; Brown et al. 2007, 2008, 2009; Pietu et al. 2007; Es- 
paillat et al. 2008). The properties of systems with gaps 
and holes are beginning to provide valuable insight into 
the physical mechanism(s) responsible for the dispersal 
of the circumstellar disk. Photophoresis, the presence of 
unknown binary companions, grain growth, dynamical 
interactions with giant planets in formation, and photo- 
evaporation have all been suggested as clearing mecha- 
nisms; there is some indication that different processes 
may dominate at different ages (see, e.g., Sicilia-Aguilar 
et al. 2010). 
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The Formation and Evolution of Planetary Systems 
(FEPS) Legacy program on the Spitzer Space Telescope 
(Meyer et al. 2006; Carpenter et al. 2008) is a spectropho- 
tometric survey of nearby sun-like stars, with masses 
from 0.7 to 1.3 Mq and ages between 3 Myr and 3 Gyr. 
These ages bracket the period of time when gas and dust 
were cleared from the primordial Solar nebula, and the 
epoch when the Solar system achieved its present con- 
figuration. For 309 objects in the FEPS sample, the 
survey includes IRAC 3.6-8.0 /xm photometry to probe 
for hot, dusty analogues to the asteroid belt in the Solar 
system, MIPS 24 and 70 /im photometry to probe dust 
in the Kuiper Belt regions, and IRS 5-40 /xm spectra to 
search for miner alogical features. The excellent cover- 
age of the infrared spectral energy distribution (SED) 
permits modeling of the temperature, size, composition, 
and an initial estimate of the spatial distribution of dust 
grains (see, e.g., Kim et al. 2005; Bouwman et al. 2008; 
Cortes et al. 2009). The IRS spectra are particularly use- 
ful for identifying systems with inner cavities or gaps in 
their dust distribution. However, there are many degen- 
eracies inherent in the derivation of spatial information 
from unresolved spectra, and the SED provides little in- 
formation about the large grains that comprise most of 
the dust mass in the system. It is therefore desirable to 
combine the information from the SED with spatially re- 
solved observations at millimeter wavelengths. Spectral 
line observations of low- level rotational emission from the 
CO molecule can provide a complementary probe of the 
molecular gas content, which is the dominant mass con- 
stituent, and can yield important clues to the evolution 
of transitional objects. 

In this paper, we analyze the SEDs and resolved mil- 
limeter continuum emission of two sources from the 
FEPS sample, RX J1842.9-3532 and RX J1852.3-3700. 
These sources were detected in the i^O^^T All- Sky Sur- 
vey and identified as young stars by Neuhauser et al. 
(2000). They have spectral type K2 and K3, respec- 
tively (Carpenter et al. 2008), and have been classified 
as classical T Tauri stars (cTTSs) based on the presence 
of strong Ha emission lines. Neither was identified as a 
multiple-star system in the speckle- interferometric obser- 
vations of Kohler et al. (2008). They are located within a 
few degrees of the CrA molecular cloud (distance 130 pc; 
Neuhauser et al. 2000), and have kinematics and secu- 
lar parallaxes consistent with the RCrA association (E. 
Mamajek, private communication). The dust grain com- 
position of the inner disks of both systems was modeled 
in detail by Bouwman et al. (2008), and Pascucci et al. 
(2007) report the detection of [Ne II] emission likely aris- 
ing from a hot disk atmosphere and calculate accretion 
rates for both systems of order lO~^M0yr~^. The es- 
timated stellar ages are ~10Myr, among the oldest in 
the 1-10 Myr range for T Tauri stars in CrA measured 
by Neuhauser et al. (2000). These sources were selected 
for study on the basis of their age, their brightness in the 
1.2mm continuum (Carpenter et al. 2005), and their ac- 
cessibility to the Submillimeter Array (SMA)^^, so that 
the spectral information from the FEPS survey could 

^^ The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics and is funded by the Smith- 
sonian Institution and the Academia Sinica. 



be combined with resolved submillimeter observations. 
They were also chosen for their proximity on the sky, 
which increases the efficiency of submillimeter observa- 
tions by allowing them to share calibrators. The high- 
quality Spitzer IRS spectra provide constraints on the 
dust composition and temperature structure on the two 
systems, and both exhibit a flux deficit in the mid-IR 
photospheric excess that points to the presence of an in- 
ner hole or gap. 

Sections 2 and 3 describe the collection of new data 
from the SMA and the Atacama Submillimeter Telescope 
Experiment (ASTE) that complement the spectra and 
broadband fluxes from the FEPS survey and the litera- 
ture. In Section 4.1 we present the tools and techniques 
that we use to model the SED and resolved millimeter- 
wavelength data, and we present the models of the dust 
disk structures in Section 4.2. In Section 4.3, we explore 
the dust disk model in the context of the constraints on 
the gas content of the two systems. We summarize our 
results and discuss their implications in Section 5. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. SMA Observations 

The SMA observations of the two sources took place 
on 2005 May 14 during a full six-hour track with six of 
the 6-meter diameter antennas operating in the compact- 
north conflguration, yielding baseline lengths between 10 
and 180 meters (8 and 140 kA). The phase stability was 
adequate for most of the track, with phase differences of 
20-30 degrees between calibrator scans, but the phases 
lost coherence during the last hour of the night. The 
weather was fair, with the 225 GHz atmospheric opacity 
increasing from 0.10 to 0.14 throughout the night. Ob- 
servations of the two FEPS sources were alternated with 
observations of the quasar J1924-292 at 15-minute inter- 
vals to calibrate the atmospheric and instrumental vari- 
ations of phase and amplitude gain. Callisto was used 
as the flux calibrator, with a calculated brightness of 
6.59 Jy; the derived flux of J1924-292 was 5.4 Jy, with 
an estimated uncertainty of < 20%. The correlator was 
configured to provide a spectral resolution of 512 chan- 
nels over the 104 MHz bandwidth in the chunk contain- 
ing the C0(2-l) line, corresponding to a velocity reso- 
lution of 0.26 km s^^. Two other 104 MHz chunks were 
observed at high resolution at frequencies corresponding 
to the 13C0(2-1) and 0^80(2-1) lines. The remainder of 
bandwidth in each 2 GHz- wide sideband was devoted to 
measuring the 230 GHz continuum, observed at a spec- 
tral resolution of 4.2 kms""'^. The data were calibrated 
using the MIR software packagC"'^^, and imaging was car- 
ried out with the MIRIAD software package. 

2.2. ASTE Observations 

Observations of both FEPS sources took place on 2008 
June 28 and 29 using the CATS345 receiver on the 10.4- 
meter ASTE dish. RX J1842.9-3532 was observed on 
both nights, while RX J1852. 3-3700 was observed only 
on the second night. The receiver was tuned to place 
the CO(3-2) rotational transition in the lower sideband 
with the HCO"*'(4-3) transition in the upper sideband. 
The high resolution spectrometer mode was used to par- 

^'^ http://cfa-www.harvard.edu/~cqi/mircook.htrnl 
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tition the 128 MHz bandwidth into 1024 channels, yield- 
ing a spectral resolution of 0.11 km/s. Position switch- 
ing was used to subtract the instrumental and sky back- 
ground. In order to ensure that the detected CO emission 
originated at the position of the star, we used an offset 
position 1.5arcmin to the east on the night of June 28 
and l.Sarcmin to the west on June 29 and averaged the 
baseline-subtracted spectra to create the final spectrum. 
The double-peaked CO(3-2) line from RX J1842.9-3532 is 
detected independently on both nights using the different 
offset positions, which allows us to localize the emission 
to within 1.5 arcmin of the star (^4 beam widths). 

The CO (3-2) and HCO+(4-3) tuning of the receiver 
resulted in a spurious 1.11 MHz sinusoidal ripple of vari- 
able amplitude and phase across the bandpass, which was 
subtracted individually from each 10-second integration 
in the following manner. The amplitude and frequency 
of the ripple were estimated by finding the peak in a 
fast fourier transform of the spectrum, and then a least 
squares fit was performed to determine the precise am- 
plitude, phase, and frequency of the sinusoidal ripple, 
plus the slope and intercept of a linear component to re- 
move the worst of the baseline features. During this pro- 
cess, the region of the spectrum containing the line was 
not included in the fit so as to avoid inadvertently sub- 
tracting it. Integrations with an abnormally large ripple 
amplitude or highly irregular baseline shape across the 
bandpass were discarded (roughly 10% of each data set). 
After subtracting the sinusoid and linear fit, a third-order 
polynomial was subtracted from each integration using 
the CLASS software package^**, fitting the 20kms^^ to 
each side of the region that appeared to contain the line. 
The 10-second integrations were then averaged together 
to produce a spectrum for each night; the spectra for each 
night were averaged and weighted according to their total 
integration time (yl/^int) to produce the final spectrum. 
Due to differing exposure times and poorer weather on 
the night of June 29, the rms noise in the spectrum of 
RX J1842.9-3532 is 24 mK, while the rms noise in the 
spectrum of RX J1852.3-3700 is 39 mK. To derive the 
absolute flux scale, we calculated the main beam effi- 
ciency using observations of the calibrator M17SW taken 
throughout the night. Assuming a peak main beam tem- 
perature in the CO(3-2) line of 85. 3K, derived on the 
CSO 10.4m telescope by Wang et al. (1994), we derive 
main beam efficiencies that varied between 0.55 and 0.63 
over the course of the two nights. 

3. RESULTS 

3.1. Millimeter Continuum 

With the SMA observations, we detect 230 GHz con- 
tinuum emission from the outer disks around both target 
stars. The contour maps in the left panels of Figure 1 
show the strong detection of emission centered on the 
star position at the middle of the field. The spatial reso- 
lution of the data is I'.'Oxl'.'T, or ^150 AU at a distance 
of 130 pc, well matched to the typical size of a circum- 
stellar disk but insufficient to provide evidence for the 
presence or absence of an inner cavity. In order to esti- 
mate the integrated flux and approximate outer disk ge- 
ometry, a Gaussian flt to the visibilities was performed 

14 http://www.iram.fr/IRAMFR/GILDAS 



using the MIRIAD task uvf it. For RX J1842.9-3532, 
the fit yields an integrated flux of 49±8mJy and indi- 
cates that the disk is only marginally resolved: the ma- 
jor and minor axes lengths of 0'.'74±0'.'32 and 0'.'44±0'.'32 
suggest that the disk is probably not viewed face-on, but 
do not place strong constraints on the inclination. We 
estimate an inclination angle of 54° based on these mea- 
surements, which is consistent with the morphology of 
the ASTE CO(3-2) spectrum described in Section 3.2 
below. The fit to the RX J1852.3-3700 visibflities yields 
a flux of 60±8 mjy and major and minor axes of length 
0'.'76±0'.'21 and 0'.'73±0'.'21, consistent with a nearly face- 
on geometry; we use these measurements to estimate an 
inclination of 16°. The inclination estimates are highly 
uncertain, but the intermediate geometry of RX J1842.9- 
3532 is supported by the line profllc in Section 3.2 and 
the more face-on geometry of RX J1852. 3-3700 is sup- 
ported by the Ha line profile modeling of Pascucci et al. 
(2007). While observations at higher resolution would be 
advantageous for constraining the detailed mass distri- 
bution, even the rudimentary estimates of disk geometry 
provided by these observations are useful for constrain- 
ing the disk properties when combined with constraints 
from the broadband SED. Simultaneous modeling of the 
SED and millimeter-wavelength visibilities is described 
in Section 4.1 below. 

3.2. C0(2-l) and C0(S-2) Line Observations 

We do not detect molecular gas emission from ei- 
ther system in the interferometric SMA observations 
of the C0(2-l), "C0(2-l), or 01^0(2-1) lines. The 
data provide a 3(7 upper limit of 0.4 Jybeam^^ in each 
0.26 km s~^ channel, with a synthesized beam size of 
l'.'5x0'.'8. Although the disks are only marginally re- 
solved in the continuum emission, there is reason to ex- 
pect that the extent of C0(2-l) emission may be sev- 
eral times larger than that of the continuum (see, e.g., 
Hughes et al. 2008b). As a result, spatial filtering may 
be a factor in the non-detection (see further discussion 
in Section 4.3). 

We do not detect any CO (3-2) emission in the ASTE 
observations of the disk around RX J1852. 3-3700, with 
an rms of 39 mK in each O.lkms"^ channel. Obser- 
vations of the disk around RX J1842. 9-3532 reveal a 
double-peaked line profile, shown as a solid black line 
in Figure 2. The integrated strength of the CO(3-2) 
line is 0.24Kknis~^ with a peak main-beam brightness 
temperature of 130mK and FWHM of 2.6kms-i. The 
double-peaked profile is consistent with material in Kep- 
lerian rotation about the star, viewed at an intermediate 
inclination angle of ^54°. We detect no emission from 
the CrA molecular cloud near the line in velocity space, 
although it is possible that absorption from the cloud in 
the vicinity of the disk might infiuence the line shape de- 
spite the sources' large distance from CrA cloud center. 
In Section 4.3 below, we investigate the relationship of 
the CO(3-2) emission to the dust properties, including 
implications for the disk geometry and gas-to-dust mass 
ratio. 

4. ANALYSIS 

In order to characterize the basic properties of the 
disks, we seek a model that can reproduce the obser- 
vational features of each system. We assemble a data 
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Fig. 1.— SMA 230 GHz maps of the continuum emission from RX J1852. 3-3700 (top) and RX J1842. 9-3532 (bottom). Each row shows three 
panels: the SMA data (left), the representative model (center), and the residuals (right). The contour levels are [2,4,6,...] x3.5mjy beam" -"^ 
(the rms noise), with solid lines indicating positive contours and dotted line indicating negative contours. The l'.'0xl'.'7 synthesized beam 
is indicated by the ellipse in the lower left corner. 



set that combines the mihinieter-wavelength properties 
of the gas and dust described above with constraints from 
the broadband SED and IRS spectrum. We use the IRS 
spectrum and SED from the FEPS database (described 
in Carpenter et ah 2008), including optical, near-IR, and 
millimeter fluxes collected from the literature (Neuhauser 
et al. 2000; Skrutskie et al. 2006; Carpenter et al. 2005). 
Figures 3 and 4 show the SED (black points) and the 
IRS spectrum (red line) in the left panel for each disk, 
alongside the SMA 230 GHz visibilities (black points) in 
the right panel. In order to improve the signal-to-noise 
ratio of the plotted data, the visibilities have been depro- 
jected (see, e.g.. Lay et al. 1997) according to the disk 
geometry inferred in Section 3.1 and averaged in bins of 
15 kA. For a mathematical description of the abscissae of 
the visibility plots, refer to Hughes et al. (2008b). 

4.1. Modeling the SED and Millimeter Visibilities 

In an effort to reproduce these observations, we gen- 
erated synthetic broadband SEDs, Spitzer IRS spectra, 
and millimeter continuum visibilities using the radiative 
transfer method and disk structure models described by 
Andrews et al. (2009). In these flared, axisymmetric 
disk structure models, the radial surface density profile 
is characterized by a similarity solution for viscous accre- 
tion disks, E oc {Rc/R)'' exp[—{R/Rc)'^~'^], where Re is 
a characteristic radius and the normalization is propor- 
tional to the disk mass (for simplicity, the radial index 
has been fixed to 7 = 1; Lynden-Bell & Pringle 1974; 
Hartmann et al. 1998). Vertically, the densities are dis- 
tributed as a Gaussian with a scale height that varies as 
a power-law with radius, H ex R^^^ , with scale height 



iJioo at a radius of 100 AU. This parametric definition 
of the vertical dust distribution is maintained to mimic 
the sedimentation of dust grains below the disk atmo- 
sphere (e.g., Dullemond & Dominik 2004b); no attempt 
is made to iterate on the density structure to force the 
dust into vertical hydrostatic equilibrium. To model the 
cleared inner disks for these transitional sources, we scale 
down the surface densities by a factor (5s inside a radius 
i?cav (Scav = (^sS; Audrcws et al. 2009). Moreover, in 
an effort to better reproduce the detailed shape and solid 
state features in the IRS spectra, we permit a small (mul- 
tiplicative) increase in the scale-height at the cavity edge 
(Sh) and adjust the dust grain properties in the inner 
disk (for details, see Andrews et al. 2010). Inner disk 
silicate opacities are taken from the Jena database (see, 
e.g., Ossenkopf & Henning 1994). 

For a given parametric disk structure, fixed input stel- 
lar information (Carpenter et al. 2008), and opacities 
(see Andrews et al. 2009), we use the two-dimensional 
Monte Carlo radiative transfer code RADMC (Dullemond 
& Dominik 2004a) to calculate an internally-consistent 
temperature structure and generate synthetic data prod- 
ucts that can be compared to the observations. How- 
ever, the parameter degeneracies introduced by the ad- 
ditional inner disk parameters and the high quality of the 
IRS spectra (requiring attention to mineralogy) make the 
minimization method described by Andrews et al. (2009) 
prohibitive. Instead, we aimed to find a representative 
model that can reproduce the basic features of the data 
by focusing on varying parameters like the cavity size 
(^cav) and surface density reduction ((5s)- These models 
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Fig. 2.— ASTE CO(3-2) spectra of the disks around RX J1852.3- 
3700 (upper) and RX J1842. 9-3532 (lower). No emission is detected 
from the RX J1852.3-3700 system. The RX J1842.9-3532 emission 
(sohd hne) displays the characteristic double-peaked profile char- 
acteristic of an inclined structure in Keplerian rotation about the 
central star. The line profile predicted by the SED- and visibility- 
based model of the dust disk structure (dotted line) compares fa- 
vorably with the observations. 

serve as initial estimates of the disk structures that can 
be refined when future observations are available (e.g., 
high angular resolution millimeter data). 

4.2. Representative Models 

Table 1 presents the parameters of representative disk 
structure models capable of reproducing the observa- 
tional data for both systems, and indicates those pa- 
rameters that were fixed by particular observational con- 
straints. The 130 pc distance to the RCrA association is 
from Neuhauser et al. (2000) and the visual extinction 
and stellar properties are drawn from the FEPS database 
(Carpenter et al. 2008), while the inclination and posi- 
tion angle are estimated from the data as described in 
Sections 3.1 and 3.2 above. The other parameters are 
defined in Section 4.1; Andrews et al. (2009, 2010) in- 
clude extensive discussion of the degeneracies between 
parameters and the ways in which the observational fea- 



tures are linked to the components of the disk structure 
model. Here we include comments on several parameters 
that are particularly relevant for reproducing the data 
described in this paper. 

Inner disk structure (Rcav o.nd Scav) — These parame- 
ters are tied primarily to the wavelength and magnitude 
of the rise in the infrared flux longward of the 10 fiiri sil- 
icate feature. It should be noted that while the density 
reduction 6cav is greater for RX J1842. 9-3532 than for 
RX J1852. 3-3700, the initial difference in surface density 
must be taken into account: because RX J1842. 9-3532 is 
almost a factor of four more compact than RX J1852.3- 
3700, the surface density throughout most of the disk, in- 
cluding within the cavity, is larger. As a result, the inner 
disk of RX J1852. 3-3700 is entirely optically thin, while 
that of RX J1842.9-3532 includes both optically thick 
and optically thin regimes. The surface density profile 
of the two models is plotted in Figure 5. As indicated 
in Section 5 below, the details of the inner disk struc- 
ture are not well constrained by these models, although 
the presence of an inner cavity of greatly reduced surface 
density is firmly indicated. 

Puffing at inner edge of outer disk (SH) — The param- 
eter (5H, introduced in Andrews et al. (2010), describes 
the extent to which the scale height at the edge of the 
cavity is puffed up, which occurs because of the frontal 
illumination of the inner edge of the outer disk and is 
primarily manifested observationally in the shape of the 
far-IR SED. While a small dH can help to account for 
the very steep mid- to far-IR jump in fiux observed in 
the RX J1852. 3-3700 IRS spectrum, no shadowing is re- 
quired to reproduce the spectrum of RX J1842. 9-3532. 

Inner disk dust properties — We have very little in- 
formation about the composition of dust grains in the 
outer disk, other than that millimeter-size grains must 
be present to account for the emission at long wave- 
lengths. However, reproducing the features in the IRS 
spectrum requires us to vary the composition of the 
hot grains in the inner disk and wall that give rise to 
these features, as in Andrews et al. (2010). The 10 /im 
silicate feature and steep rise in fiux near 20 fim from 
the disk around RX J1852. 3-3700 are reproduced well 
by an inner disk and cavity edge populated entirely 
by small (~0.1/im) amorphous silicate grains. For the 
disk around RX J1842. 9-3532, the strength and posi- 
tion in wavelength of the silicate feature are well re- 
produced by a mixture of small (amax=l/^ni) and large 
(amax=l mni) ISM-composition grains (^^80%) and crys- 
talline and amorphous silicates (~20%) in the inner disk. 
This combination of grain compositions is by no means 
a unique solution to the problem of fitting the mid-IR 
spectrum, but merely demonstrates that a mixture of 
different grain properties is helpful in accounting for the 
observed spectral features. A more detailed mineralogi- 
cal analysis of these systems can be found in Bouwman 
et al. (2008). The estimated mass of dust in the inner 
disk is 10-1° M0 for RX J1842.9-3532 and 10"^ M© for 
RX J1852. 3-3700. 

The model SED and millimeter visibilities for the 
structural parameters in Table 1 are shown by the green 
lines in Figures 3 and 4. They reproduce the basic fea- 
tures of all of the available dust disk diagnostics, in- 
cluding the broadband SED, the IRS spectrum, and the 
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Fig. 3. — Spectral energy distribution (left) and the real and imaginary components of the deprojected SMA 230 GHz visibilities (right) 
for RX J1842. 9-3532. The broad-band SED (black points) and IRS spectrum (red line) are well reproduced by the best-fit RADMC disk 
structure model (green line). The model stellar photosphere (dashed blue line) is plotted for comparison. The units of the ordinate are 
defined so that L^ = incPuFu in units of Lq . For a mathematical definition of the abscissa, refer to Hughes et al. (2008b); the deprojection 
is carried out as in Lay et al. (1997). 



RX J1852. 3-3700 

Model 

IRS Spectrum 

Photosphere - 




10 100 1000 

A (yu,m) 

Fig. 4.— Same as Fig. 3 above, but for RX J1852. 3-3700. 



80 
60 
40 
20 


20 


-20 




H — I — I — I — I — I — I — I — I — H 



<• o I 



20 40 60 80 100 120 



millimeter-wavelength visibilities. In the discussion be- 
low, we focus on the most robustly-constrained model 
parameters, including the extent and surface density re- 
duction of the inner cavity and the size and dust mass of 
the outer disk. 

4.3. Constraints on Molecular Gas Content 

Here we compare the predictions of the dust disk model 
with the constraints on the CO emission described in 
Section 3.2. For simplicity, we assume that gas and 
dust are well-mixed, with identical spatial and tempera- 
ture distributions. We use the standard assumptions of 
a gas-to-dust mass ratio of 100:1 and a CO abundance 
of lO"** relative to H2. To account for the depletion of 
gas-phase CO onto dust grains at low temperatures, we 
decrease the CO abundance by an additional factor of 
lO""* in regions with temperatures <20K. As in Andrews 
et al. (2009), we then use the Monte Carlo molecular line 
radiative transfer code RATRAN (Hogerheijde & van der 



Tak 2000) to calculate the level populations and predict 
the sky-projected intensity of CO arising from each sys- 
tem, given the underlying structure of the representative 
models derived in Section 4.2. We use the MIRIAD task 
convol to convolve the resulting intensity distributions 
with the 21'.'1 beam of the 10.4 m ASTE telescope, since 
the ASTE spectra provide the most stringent limits on 
the CO emission from the systems. 

With the standard assumptions of gas-to-dust and CO- 
to-H2 ratios, the model projection of CO (3-2) line flux 
strongly overpredicts both the upper limit for the disk 
around RX J1852. 3-3700 and the detection of CO(3-2) 
emission from the disk around RX J1842. 9-3532 by a 
factor of 4.8. In order for the model to successfully re- 
produce the weak emission from RX J1842. 9-3532, the 
volume density of CO must be decreased to 8±3% of its 
initial value. A model spectrum for this case is given 
by the dotted line in Figure 2, and compares well with 
the strength and width of the observed CO spectrum. 
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TABLE 1 

Star and Estimated Disk Parameters 



Parameter 


RX J1842.9-3532 


RX J1852. 3-3700 


Fixed 


Distance (pc) 


130 


130 


Ay (magnitudes) 


1.06 


0.97 


Spectral Type 


K2 


K3 


Logg (logcms-2) 


4.27 


4.31 


Toff (K) 


4645 


4854 


in 


54 


16 


P.A. {°) 


32 


-56 


7 


1.0 


1.0 


Varied 


Rc (AU) 


50 


180 


md (M@r 


0.010 


0.016 


^ 


0.2 


0.2 


Hloo (AU) 


4.8 


6.3 


Rcav (AU) 


5 


16 


<5cav 


9 X 10-6 


3 X 10-"^ 


<5H 


1 


1.4 



''Total mass in gas and dust, assuming a gas-to-dust mass ratio 
of 100 



While the line peaks appear narrower than the model, 
the noise in the line is too large to merit modeling the 
profile in detail; it is also possible that contamination 
from remnant molecular cloud material could contribute 
to the narrowing of the peaks. 

It is worthwhile to consider whether the apparently 
faint CO emission from these disks could be produced by 
overestimated model gas temperatures, rather than low 
CO abundances. In our models, we have assumed that 
the gas and dust are well-mixed. However, the infrared 
emission from these sources is too faint to accomodate a 
dust distribution in vertical hydrostatic equilibrium: the 
dust has apparently settled toward the disk midplane. 
The gas, on the other hand, should still feel a verti- 
cal pressure gradient, and therefore have a larger scale- 
height than the dust (e.g., Dullemond & Dominik 2004b; 
D'Alessio et al. 2006). As a result the gas in our model, 
which is assumed to follow the spatial distribution of the 
dust, is at somewhat lower temperatures than it would be 
in a thicker disk. This difference in temperature due to 



the physical location of the gas and dust is apparently ex- 
acerbated by exposure to high-energy (FUV, X-ray) radi- 
ation from the central star (Qi et al. 2006). Simulations 
show that these heating processes typically raise the gas 
temperature over that of the dust in the upper disk layers 
from which the CO emission arises (e.g., Kamp & Dulle- 
mond 2004; Woitke et al. 2009). If anything, our model 
assumptions should then bias the results toward underes- 
timating gas temperatures in the disk atmospheres, and 
therefore under-predict the amount of CO emission. 

If instead we take the well-mixed models at face value 
(i.e., assuming identical gas and dust temperatures), the 
model gas temperatures can only be modified by adjust- 
ing the dust structure parameters. The model tempera- 
tures in the outer disk where the CO emission is gener- 
ated are primarily set by the size of the disk (i.e., distance 
from the star) , and the amount of flaring. To reduce the 
gas temperatures by the factor of ^^5 required to repro- 
duce the observations for standard CO/II2 and gas-to- 
dust conversions, we would need to increase the emitting 
area {Rc) by roughly an order of magnitude or decrease 
the scale-height {H) to an extreme, flat-disk level: these 
possibilities are inconsistent with the millimeter visibili- 
ties and infrared SED, respectively. Given these difficul- 
ties, a more straightforward interpretation of the faint 
line emission in these disks is that the CO abundances 
are substantially lower than would be expected for the 
standard gas-to-dust mass ratio of 100 and/or CO abun- 
dance relative to II2 of 10~* (see also Qi et al. 2004; 
Hughes et al. 2009). 

In the absence of measurements of the H2 content of 
the outer disk, it is not possible to determine whether 
the lower CO content results from a reduction in gas-to- 
dust mass ratio or abundance of CO relative to H2, but in 
either case it marks a significant departure from standard 
assumptions. Because the rms in the ASTE spectrum 
of RX J1852.3-3700 is larger than that of RX J1842.9- 
3532, the upper limit on the CO content of the disk is 
similarly ~8% of the initial value, assuming a standard 
gas-to-dust ratio and CO abundance. The model with 
reduced CO content relative to standard assumptions is 
also consistent with the limits on C0(2-l) emission from 
the SMA. We use RATRAN to generate a sky-projected 
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C0(2-l) emission map, which is then sampled with the 
fourier components of the SMA data using the MIRIAD 
task uvmodel to account for spatial filtering effects. The 
model with standard CO abundance should be detected 
by the SMA observations, whereas the model with 8% 
CO content is consistent with the upper limits on the 
C0(2-l) emission from both systems. 

If the gas disk were truncated relative to the dust disk, 
this could contribute to the low CO content. It is un- 
likely, however, that both systems would undergo trun- 
cation - especially given the dearth of companions within 
6" (Kohler et al. 2008) - and the truncation would have 
to be severe in order to account for an order of magni- 
tude reduction in CO content. It should also be noted 
that the conclusion of reduced CO content is largely in- 
dependent of the model parameters describing the inner 
disk and the transitional region between inner and outer 
disk. The CO(3-2) emission arises only from the cold 
outer disk, and the gas-to-dust ratio is derived only for 
this region, since we have no information on the gas con- 
tent of the inner disks of these systems. The extent of 
the outer disk and its total dust mass are derived from 
two observational parameters: the millimeter flux and 
size scale indicated by the resolved visibilities. While the 
vertical structure and inner disk properties can affect the 
temperature of the outer disk and therefore the magni- 
tude of the CO (3-2) emission, these effects are secondary 
to the basic midplane temperature structure determined 
by the radial scale of the dust disk. The reduced CO con- 
tent is therefore robust to variations in the details of the 
inner disk structure, since variations in inner disk proper- 
ties will have only second-order effects on the gross outer 
disk properties from which this conclusion is derived. 

Both systems therefore appear to have undergone a re- 
duction in molecular gas content relative to the standard 
assumptions for primordial disks. Given their age and 
transitional SEDs, this may indicate that gas dispersal 
is underway simultaneously with dust clearing from the 
inner disk. 

5. DISCUSSION 

We have generated models that can reproduce simul- 
taneously the basic observational diagnostics of the gas 
and dust disks around RX J1842.9-3532 and RX J1852.3- 
3700, including their broadband SEDs, IRS spectra, re- 
solved millimeter- wavelength visibilities, and CO (3-2) 
spectra. As indicated by the mid-IR flux deficit, both 
systems are transitional, with an inner cavity of signifi- 
cantly decreased dust optical depth. 

The disk around RX J1842.9-3532 also exhibits a sub- 
stantial near-IR excess over the stellar photosphere. It 
shares this feature with the sample of objects labeled 
gapped, or "pre-transitional" by Espaillat et al. (2007). 
They model such systems with an optically thin inner 
disk bracketed by an optically thick ring close to the star 
and the optically thick outer disk at large radii. Similar 
models for the LkCa 15 system, refined with the addition 
of radiative transfer through the inner disk, are described 
in Mulders et al. (2010). IseUa et al. (2009) modeled 
the near-IR excess and mid-IR deficit in the LkCa 15 
system using a density distribution that increases with 
distance from the star, but includes a puffed-up inner 
rim at the dust disk edge. In our study, the inner disk 
model retains the continuous surface density profile of the 



outer disk (decreasing with distance from the star), sup- 
pressed by the factor Scav (see Figure 5), with no change 
in scale height at the inner edge of the inner disk. Due 
to the relatively small scale heights in the inner disk, 
we can approximate the optical depth to starlight as the 
product of surface density and 1 /xm opacity, ErKi^hi. 
In this approximation, the cavity is optically thick be- 
tween 0.01 and ^0.2 AU but optically thin between ~0.2 
and 5 AU, comparable to the models described in Espail- 
lat et al. (2007). These results suggest that transition 
disks with near-IR excess are not necessarily "gapped" 
in terms of their surface density or discontinuous in terms 
of their scale height, since we demonstrate that the inner 
disk can be modeled using a single, continuous surface 
density function for the disk cavity that contains just 
enough mass to have both optically thick and optically 
thin regimes. Effectively, this indicates that we can place 
no constraint on the contrast in surface density or scale 
height between the "gap" and the optically thick ring 
near the star based on the morphology of the IRS spec- 
trum. This is refiected by the success of several very 
different models of inner disk structure (Espaillat et al. 
2007; IscUa et al. 2009; Mulders et al. 2010, this work) in 
reproducing the characteristic mid-IR deficit surrounded 
by near- and far-IR excesses. 

The model parameters for the two disks discussed in 
this paper are consistent with those of the nine disks 
in Ophiuchus that were studied using this method by 
Andrews et al. (2009). This is perhaps unsurprising, 
since these targets were similarly selected on the basis 
of their large submillimeter fluxes. RX J1842. 9-3532 and 
RX J1852. 3-3700 have slightly lower masses due to miss- 
ing material in the otherwise dense disk center, as for 
the transitional systems in the high-resolution Ophiuchus 
sample. Yet as with the other transition disks in Andrews 
et al. (2009) they are still on the high end of the distri- 
bution of masses of Taurus and Ophiuchus disks in the 
sample of Andrews & Williams (2005, 2007). The tran- 
sition disk frequency for very low-mass stars in the CrA 
region has been examined by Sicilia-Aguilar et al. (2008); 
Ercolano et al. (2009), with estimates ranging from 15- 
50%. On the whole, the disks around RX J1842.9-3532 
and RX J1852. 3-3700 appear comparable to transition 
disks in other low-mass star forming regions. 

While relatively little is known about the gas evolu- 
tion of circumstellar disks, it is somewhat surprising that 
such massive dusty disks should have such low CO con- 
tent. It is interesting to note that other transitional 
systems around K stars have also been found to have 
low CO content through similar modeling approaches 
to that presented here: TW Hya (Qi et al. 2004) and 
CM Aur (Hughes et al. 2009) have derived CO abun- 
dances of 10^"'' and 10~^, respectively (see also Kastner 
et al. 1997). The low gas content of TW Hya has re- 
cently been corroborated by Herschel observations that 
appear to indicate a gas-to-dust mass ratio of ~2.6-26, 
substantially lower than the standard value of 100 de- 
rived for the ISM (Thi et al. 2010). In the absence of 
knowledge about the H2 content of the system, this is 
consistent with the Qi et al. (2004) measurement of low 
CO abundance, which assumed a gas-to-dust mass ratio 
of 100. To be clear, the inferred CO content of these 
systems is judged to be "low" in comparison with the 
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standard assumptions of a gas-to-dust mass ratio of 100 
and a CO abundance relative to H2 of 10^''. A reduc- 
tion in either total molecular gas content (i.e., gas dis- 
persal) or CO relative to H2 (e.g., molecular depletion) 
could be responsible, or the standard assumptions could 
be incorrect. While apparently consistent with ISM con- 
ditions, there are almost no observational constraints on 
either gas-to-dust or CO-to-H2 ratios in the outer regions 
of primordial (non-transitional) protoplanetary disks. It 
is particularly important to understand depletion, since 
molecular depletion onto dust grains has been suggested 
as a possible cause of low molecular gas content in cir- 
cumstellar disks (see, e.g., Dutrey et al. 1997). Depletion 
is sensitive to the details of the temperature and density 
structure of the disk, which may be difficult to determine 
at low spatial resolution. The high resolution survey of 
Andrews et al. (2009) finds CO abundances consistent 
with standard assumptions to within a factor of a few for 
the three Ophiuchus disks exhibiting strong CO emission, 
although they note that most of the sources in the sam- 
ple are substantially contaminated by the local molecular 
cloud within which the stars are embedded. Comparable 
modeling efforts for solar-type systems less embedded in 
their natal clouds arc therefore highly desirable for deter- 
mining the implications of the low CO content of these 
transitional systems. 

Assuming that the standard assumptions are correct 
for protoplanetary disks, it is interesting to consider the 
implications of low CO content for the physical processes 
clearing the cavities in transitional systems. Among the 
several mechanisms proposed to clear the central regions 
of disks, the most popularly invoked are gravitational in- 
teraction with a giant planet in formation (e.g., Lin & Pa- 
paloizou 1986; Brydcn et al. 1999) and photoevaporation 
by energetic radiation from the star (e.g., Clarke et al. 
2001; Alexander et al. 2006). Many of the dust prop- 
erties of transitional systems, including the disk masses 
inferred from dust flux alone and the sharp transition 
between the inner and outer disk from SED modeling, 
point towards planet formation as the most likely pro- 
cess. The mass accretion rates of the observed systems, 
particularly in combination with the disk masses from 
dust, are also generally more consistent with planet for- 
mation as the dominant physical mechanism (see, e.g., 
discussions in Najita et al. 2007; Alexander & Armitage 
2007). RX J1852.3-3700 in particular occupies an in- 
teresting and underpopulated region of the M-M^, plane 
investigated by Najita et al. (2007), due to its extremely 
low accretion rate measurement. The identification of 
gapped or "pre-transitional" systems like RX J1842.9- 
3532 has also been interpreted as being more consistent 
with planet formation than photoevaporative clearing 
(e.g., Espaillat et al. 2007). It is worth noting, however, 
that a reduced gas mass for transitional systems relative 
to protoplanetary disks would present some difficulty for 
the planet formation scenario, since there is no reason a 
priori to expect the formation of giant planets to alter 
the molecular gas content of the outer disk. The masses 
in Table 1 may be misleading, since they represent the 
total mass in gas and dust of the RADMC model assuming 
a standard gas-to-dust mass ratio of 100, without taking 
into account the evidence for low CO content described 
in Section 4.3: if the low CO content is a result of reduced 



gas-to-dust mass ratio, the disk masses may be reduced 
by an order of magnitude or more, placing them squarely 
within the region of parameter space preferred by photo- 
evaporative clearing models. If this were generally true 
of transition disks (including TW Hya and CM Aur), 
it could also provide an alternative explanation for the 
trend in mass accretion rate as a function of disk mass 
noted by Najita et al. (2007), since an order of magni- 
tude difference in total mass between the protoplanetary 
and transition disks would bring the two populations into 
agreement. 

Photoevaporation of disk material by energetic stel- 
lar radiation can similarly account for many, but not 
all, features of these systems. Its most obvious suc- 
cess is in explaining the low CO content of the outer 
disk, which cannot easily be accomplished within the 
framework of planet formation. One complicating fac- 
tor, however, is the measured accretion rate of mate- 
rial onto the stars based on modeling of the Ha pro- 
files by Pascucci et al. (2007). They calculate an ac- 
cretion rate of 1 x lO^^Moyr^i for RX J1842. 9-3532 
and 5 x lO^^^Mgyr-i for RX J1852.3-3700. These are 
roughly an order of magnitude below the average for 
IMyr-old stars in Taurus (GuUbring et al. 1998; Calvet 
et al. 2004), which is consistent with the trend for tran- 
sitional systems in Taurus noted by Najita et al. (2007). 
RX J1852. 3-3700 in particular occupies an underpopu- 
lated region of the M-M* plane investigated by Najita 
et al. (2007), due to its extremely low accretion rate mea- 
surement, which could indicate the presence of a very 
massive planet in the disk. This relatively low but mea- 
surable accretion rate is inconsistent with the original 
predictions of most photoevaporation models (see, e.g., 
Alexander & Armitage 2007). Recent work by Owen 
et al. (2010) suggests that accretion rates of this mag- 
nitude may in fact be consistent with the predictions of 
radiation-hydrodynamic models that incorporate both x- 
ray and EUV photoevaporation; however, they assume a 
stellar x-ray luminosity far greater than is typically ob- 
served, and EUV luminosities of young stars are not yet 
well determined. Another potential problem is that the 
10^ yr timescales predicted for photoevaporative clearing 
is only marginally consistent with a ^10% transition disk 
fraction at ages of 10^-10^ yr. Certainly more observa- 
tions are needed to understand the physical mechanisms 
responsible for clearing the inner cavities of transition 
disks, including observations across a larger sample and 
a range of ages, to determine the molecular gas content 
of transitional relative to primordial disks and constrain 
the processes responsible for their dissipation. 

At least two well-studied systems at ages of ~10Myr 
have been identified as candidate transition disks under- 
going photoevaporation: the M star HD 100453 (Collins 
et al. 2009) and the A star 49 Ceti (Hughes et al. 2008a). 
The former exhibits a strong IR excess indicative of an 
optically thick outer disk, but with no evidence of accre- 
tion and a maximum gas-to-dust ratio of 4:1 in the outer 
disk. The latter exhibits dust properties similar to a 
debris disk, yet retains an extended optically thin molec- 
ular gas disk with an inner hole. The differences between 
these systems and the CrA transition disks in this paper 
are striking, and are perhaps indicative of the range of 
evolutionary paths over which the transition from proto- 
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planetary to debris disk may occur. 

6. SUMMARY AND CONCLUSIONS 

The lOMyr old systems RX J1842.9-3532 and 
RX J1852. 3-3700 present a rare opportunity to study the 
late-stage evolution of circumstellar disks. We assemble a 
multiwavelength data set including the broadband SED, 
Spitzer IRS spectrum, and generate models of the disk 
structure for each system that can reproduce the fea- 
tures of the dust continuum emission. We infer the pres- 
ence of an optically thin 16 AU cavity in the disk around 
RX J1852. 3-3700. The near-IR photospheric excess and 
mid-IR deficit in the SED of RX J1842.9-3532 allow us 
to infer the presence of an optically thin 4AU gap in 
the disk around RX J1842. 9-3532, bracketed by an opti- 
cally thick ring near the central star (^0.01-0.2 AU) and 
an optically thick outer disk. This configuration marks 
it as a member of the recently identified class of "pre- 
transitional" objects (Espaillat et al. 2007). Comparison 
of the dust disk models with spatially unresolved C0(3- 
2) spectra from ASTE allow us to infer a molecular gas 
content roughly an order of magnitude less than that ex- 
pected for standard assumptions about the gas-to-dust 
mass ratio and CO abundance relative to H2. If the stan- 
dard assumptions are correct for primordial disks, then 
this low CO content is suggestive that evolution of the 
gas disk may be underway simultaneously with the dust 
disk dispersal. 

The transitional systems in CrA observed in this paper 
therefore mark interesting test cases for distinguishing 
between proposed mechanisms for gas and dust clearing 
at late ages. Follow-up of these objects with instruments 
at various wavelengths can help to fill in our picture of 
the properties of the inner and outer disks. For example, 
the far superior spatial resolution of the Atacama Large 



Millimeter Array (ALMA) will permit vastly improved 
modeling of the structure of the extended gas and dust 
disk, as well as providing direct access to the conditions 
within the cavity, removing ambiguity about surface den- 
sities and scale heights in the inner disk. The sensitivity 
to spectral line emission provided by ALMA and Her- 
schel will yield insight into the gas mass and chemistry 
and therefore the origin of the reduced CO content of 
the outer disk. In the meantime, observations of rovi- 
brational lines can aid in determining the gas content 
of the warm inner disk, which will aid in distinguishing 
between proposed clearing mechanisms. Scattered light 
images would also be useful for constraining the vertical 
structure of the disks and reducing degeneracies in these 
initial models. The suite of instruments currently com- 
ing online is poised to revolutionize our ability to char- 
acterize the physics of individual disks in the compelling 
transitional stage of evolution. 
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